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Abstract 

Final state radiation in the process e + e~ — ► ir + ir~ is considered for the cuts used in 
the analysis of KLOE data at large angles. By means of a Monte Carlo event generator 
FEVA, effects of non-pointlike behaviour of pions are estimated in the framework of 
Resonance Perturbation Theory. An additional complication related with the <j) meson 
intermediate state is taken into account and the corresponding contributions (the direct 
decay (j> — > 7r + 7r~7 and the double resonance decay <fi — ► /) ± 7r T — > 7r + 7r~7) are added to 
FEVA. A method to test effects of non-pointlike behaviour of pions in a model-independent 
way is proposed. 



PACS: 13.25. Jx, 12.39.Fe, 13.40.Gp 



1 Introduction 

The ongoing experiments on precise measurements of the cross section of e + e~ annihilation 
into hadrons aim to a precision at the 0.5 — 1.0% level [HOE]. Such an accuracy is crucial for 
various tests of the Standard Model [I], e.g., by confronting the experimentally measured value 
of [5], the muon anomalous magnetic moment, to the theoretical prediction. The accuracy of 
the theoretical calculation of is currently limited by the hadronic contribution, a^ ad . This 
contribution cannot be reliably calculated in the framework of perturbative QCD, because low- 
energy region dominates. Fortunately, its leading order part, be estimated from 
the dispersion relation using the experimental cross sections of e + e~ annihilation as an input 
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where <Jh{q 2 ) is the total hadronic cross section, cr^q 2 ) is the total cross section of the process 
e + e" —>■ 7 — > and the function K(q 2 ) is a smooth function that increases from 0.63 at the 

threshold (s = 4m 2 ) upto 1 for q 2 — > oo. The quantity g 2 is the total four momentum squared 
of the final hadrons. 

This behaviour of the integrand results in that the biggest contribution (about 70%) to the 
leading order hadronic part of the muon anomalous magnetic moment, o^ ad,LO \ comes from the 
energy region below or about 1 GeV. Due to the presence of the p-meson the main contribution 
is related with the 7r + 7r~ final state. 

Experimentally, the energy region from threshold to the collider beam energy is explored 
at the ^-factory DA$NE (Frascati, s = AE 2 = mj) |S1 Ej and B-factory, PEP-2 (SLAC, s = 
m T(45)) an d KEKB (KEK, s = Ww 4 m) [HI E] using the method of radiative return [TOl [TT1 [12] . 
This method relies on a factorization of the radiative cross section in the product of the hadronic 
cross section and a radiation function H(q 2 , 9 max , 9 min ) known from Quantum Electrodynamics 
(QED) [12j [131 El- The large luminosity of the $ and S-factories allows to compensate the 
a 

additional factor — caused by the hard photon emission. 

For a two-pion final state considered here it means that the radiative cross section a 71 ^ 
corresponding to the process 

e + (pi) + e-(p2) -> 7T+(p + ) + 7T(p_) + 7 (fc), (2) 

can be written as 

da™ 1 

q 2 ^ r = v 7T7r (q 2 )H(q 2 ,6 max ,6 min ), q = P+ +p_, (3) 

where the hadronic cross section a n7T is taken at a reduced CM energy. This factorization is 
valid only for photon radiation from the initial leptons (initial state radiation (ISR)). This is not 
possible for final state (FS) radiation (FSR) which is an irreducible background in radiative 
return measurements of the hadronic cross section [21 [15]. Indeed, the FSR cross section 
calculation has an additional complication compared to ISR case. In principle, RC caused by 
initial state radiation, i.e. the function H(q 2 ,9 max ,9 min ), can be calculated in QED, although 
the accuracy is technically limited. Instead, for the calculation of the FSR cross section the 
situation is different since its evaluation relies on models describing pion-photon interaction. 
Usually the combined sQED*VMD model is assumed as a model to calculate FS Bremsstrahlung 
process [T2l [T7] . In this case the pions are treated as point-like particles (the sQED model) and 
the total FSR amplitude is multiplied by the pion form factor, that is estimated in the VMD 
model. Unfortunately, the sQED*VMD model is an approximation that is valid for relatively 
soft photons and it can fail for high energy photons, i.e near the tt + tt~ threshold. In this energy 
region the contributions to FSR, beyond the sQED*VMD model, can become important. As it 
was shown in Ref. [18J the Resonance Perturbation Theory (RPT) is an appropriate model to 
describe photon-meson interactions in the energy region about 1 GeV and we will apply this 
model to estimate the Bremmstrahlung FS contributions beyond of sQED. 

At thet $ factory DA$NE there is an additional complication related with the possible 
intermediate meson state and the corresponding contributions should be included in the 
Monte Carlo event generator. 

1 In fact the process of FSR cannot be excluded from the analysis. It can be suppressed choosing the small 
angle kinematics (0 7 < 9 max -C 1) but for the large angle analysis this contribution becomes upto 40% of ISR 
and should be estimated very carefully. (For some advantages of large angle analysis compared to small angle 
one see [US].) 
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In this paper we present the results obtained by the Monte Carlo event generator FEVA 
that simulates the process ([2]) for the DA$NE accelerator setup. Our computer code FEVA 
was inspired by MC EVA [12]. The previous version of FEVA was described in [19], where the 
Bremsstrahlung process (both in the framework of RPT and sQED) and the direct decay 
(only /o parameterization) were considered. The current version of FEVA includes in addition 
the double resonance contribution and a more sophisticated parameterization for the direct 
decay. 

This paper is organized as follows. In Section [2] we give a general description of FSR process 
and present the FSR models that have been already included in our program FEVA. In Section [3] 
the numerical results for the KLOE large angle analysis are presented. Since most of the effects 
arising in the FSR are model-dependent, we conclude this paper by suggesting a way to test 
possible effects beyond sQED, in a model-independent way (see Section H]). A conclusion is 
given in Section [5j 



2 Final state radiation models 

The cross section of the FSR process can be written as 

daF = _}__ [ 5 *(q _ p p _ _ fc ) d 3 P + d 3 P-d 3 k {FSR) 2 

2s(2vr) 5 J w y+ y ! 8E + E_u 1 1 K J 

where Q = p\ + p 2 , s — Q 2 , 

m (fsr) = ^M^u(- Pl ) % u(p 2 )et, (5) 
and where the FS tensor M^ u describes the transition 

7 *(g)^ 7r +(p + ) 7 r-(p_)7(A;). 

It is convenient to parametrize the FS tensor in terms of three gauge invariant tensors (see [2U] 
and Ref. [23,24] in it): 

M^(Q, k, I) = -ze 2 Mjr(Q, k, I) = -ie\r^h + Tf/a + rffs), (6) 

= k^-g^k-Q, l=p + -p_, 
t% v = k ■ l(FQ u - g^k ■ I) + r(k"k -l-Fk- Q), 

= Q 2 (g^k ■ I - m v ) + Q»{l v k ■ Q - Q v k ■ I). 

We would like to emphasize that this expansion is totally model- independent. The model 
dependence is related only with an explicit form of the scalar functions /j (we will call them 
structure functions). 

Here is the list of the FSR processes included in FEVA MC: 

e + + e — > 7r + + 7r + 7 Bremsstrahlung process (7) 

e + + e~ — ► <fi — » (/ ; fo + cr)7 — > tt + + n~ + 7 (j) direct decay (8) 

e + + e _ — > (ft ~^ p ± vr =F — > n + + 7r~ + 7 VMD contribution (9) 

e + + e~ -> p ± 7r T -> tt + + + 7. (10) 

In the next sections we present the models describing these processes. The presence of (jSJ) and 
([H]) is due to the energy at which KLOE is running(s = mV). 
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2.1 Bremsstrahlung process 

As it was mentioned in Introduction, the sQED*VMD model is an approximation to describe 
soft photon radiation by pions. To estimate the contributions beyond the sQED*VMD model 
we applied RPT. The model is based on Chiral Perturbation Theory (%PT) with the explicit 
inclusion of the vector and axial-vector mesons, po(770) and ai(1260). Whereas gives 
correct predictions on the pion form factor at very low energy, RPT is the appropriate framework 
to describe the pion form factor at intermediate energies (E ~ m p ) [18] and satisfies the QCD 
high energy behaviour. 

Using the result of Ref. pO] we write the contribution to the functions fa (see Eq.flE])) caused 
by the Bremsstrahlung FS process as 

fi = f! QED +A f r T , (ii) 

rsQED 2fc ■ QF n (Q 2 ) rsQED _ —2F n (Q 2 ) /,„n 

h ~ (k ■ QY - (k ■ 2 ' 2 ~ (k ■ Q) 2 - (k ■ /) 2 ' 1 ] 

f.f ED = 0, (13) 



where 



Af n PT = F 2 -2F v G v fl + 



F 



2 



f 2 V m p m 2 p -s- im p T p (s) 

(k-l) 2 (s + k-Q)[4m 2 a - (s + l 2 + 2k-Q)} 



f 2 m 2 a [ D(l)D(-l) 8D(l)D(-l) 



(14) 



AfRPT _ F 2 A 4m 2 a -(s + l 2 + 2k.Q) 

Af2 ~ f*ml 8D(l)D(-l) ' l ~ P+ P ~> (15) 
= 73!«30' D { ^ m 2 a - { s + f + 2k Q + , km (16) 



For notations and details of the calculation we refer a reader to [20]. Fy, Gy and Fa are 
parameters of the model. According to the RPT model the pion form factor, that includes 
p — uj mixing, can be written as: 

FM 2 ) = 1 + ^ W) |^l " ^B^q 2 ) I , (17) 

where 



m 2 — q 2 — im r T r (q 2 ) ' 

q 2 is the virtuality of the photon, f n = 92.4 MeV and the parameter Il pcJ describes p-uo mixing. 
An energy-dependent width is considered for the p meson: 



rp(^) = W^(^-^j ' -6(^-4^), (19) 




2 In that paper it was shown that the coupling constants of the effective chiral lagrangian at the order p 4 are 
essentially saturated by meson resonance exchange. 
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while for the w-meson a constant width is used, = 8.68 MeV, and m w = 782.7 MeV. We 
assume that the parameter H pLU , that determines p-u mixing, is a constant and is related to the 
branching fraction Br{uj — ► 7r + 7r~): 

Br{u - tt+O = -ffldl (20) 

The value of i*V and Gy as well as the mass of the p meson (m p ) and the parameter of the p-u 
mixing Tl pu) were estimated by the fit of Novosibirsk CMD-2 data for the pion form factor [TJ : 

m p = 774.97 ± 1.4 MeV, n pw = -2774 ± 291.2 MeV 2 , 
T p = 145.21 ± 2.6 MeV, F v = 154.22 ± 0.5 MeV. 

Then Gy = 64.6 ± 0.3 MeV and Br(u -> tt+tt-) = (0.96 ± 0.19)%. 

For the a\ meson we take m a = 1.23 GeV and Fa = 0.122 GeV corresponding to the mean 
value of the experimental decay width F(ax — > nj) = 640 ± 246 keV [2~T] . 

We would like to mention here that the contribution of any model describing the Bremsstrahlun^ 
FS process can be conveniently rewritten as in Eq. (1111) and in the soft photon limit the results 
should coincide with the sQED*VMD model prediction. 

2.2 (f) direct decay 

For DA$NE energy (s = m 2 ) there are contributions to the final state 7r + 7r~7 related with the 
intermediate <fi meson state. In this section we consider the direct rare decay <fi — > 7r + 7r - 7. 

The direct decay is assumed to proceed through the intermediate scalar meson state 
(either f or f Q + a): <fi — > (f ; f + a) r y — > 7r + 7r~7, and its mechanism is described by a single 
form factor f^Q 2 ). As it was shown in [22], [23], this process affects the form factor /i of Eq.flH]): 

jscai _ 9^nf(j>{Q ) /2j\ 
1 s - ml + im^T^ ' 

First, we consider the case of the f intermediate state. To estimate this contribution 
we apply the Achasov four-quark model described in [21]: the <fi — > / 7 decay amplitude is 
generated dynamically through the loop of charged kaons. The form factor fs reads: 



f K+K- (0 2^ = 94>K+K- g fo n+n-gf K+K-e t5B ^ ( ^ Q_\ (OO) 

" ' " 2**m* K ( m % -Qi + ReU f0 (m%) ~ n /o (Q 2 )) Wk < ' ' 



where /(.,.) is a function known in analytic form [221 125] and 5b(Q 2 ) — by^Q 2 — Am 2 , b = 
75°/GeV. The term ReUf (mj o ) — Uf (Q 2 ) takes into account the finite width corrections to 
the f propagator [24]. A fit to the KLOE data <p — > n ^- ^^ gives the following values of the 
parameters [26J: 

m fa = 0.962 GeV, 9 f° K+K ~ = 1.29 GeV 2 , 9f f +K ~ = 3.22. (23) 

47T qj + _ 



! r(/ - n+ir-) = |r(/ - nw) 
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In a refined version of this model which includes the a meson in the intermediate state 
the form factor can be written as 

f K+K- (n 2\ _ 94>K+K~ A(&kAQ 2 )+&kk(Q 2 ))t f m l Q 2 \\^ n ^ r~ X n , 

h W ) - 9 o m 2 e 

m K \ m K m KJ RRI 

where Grr* is the matrix of inverse propagators [27] . Such an extension of the model improves 
the description of the data at low Q 2 (see Fig. 1, left) and gives the following value of the 
model parameters [28] 

2 2 

m fo = 0.977 GeV, 9 ^ K+K ~ = 1.12 GeV 2 , 9f f +K ~ = 6.9, (24) 

4tt g 2 + _ 



,2 „2 
n C 



0.462 GeV, 9 ° K+K ~ = 0.024 GeV 2 , 9a ? +K = 0.052. (25) 
An 9 + - 



2.3 VMD contribution 

Another contribution producing the intermediate meson state is the double resonance con- 
tribution (Q. In this case the off-shell meson decays to p ± and tt t followed by p — > 7r7. The 
explicit value for the functions fY MD for this decay can be found in Ref. [23]. To correspond 
to the KLOE analysis [28] we added also the additional phase between VMD and direct 
contributions, the factor Y\X MD and the phase of the cu-0 meson mixing (3^: 

jVMD = __L_^_l + 3 x + a)( ^ i)+ ^ (x2)) 

+-(x!- x 2 ){g{x 1 ) - g{x 2 )) S j, (26) 

jVMD = __±_ (g{xi)+g{x2)): 



inas 2 



where 



pn<t> nP m 2 e i Pp P i ^4> „2yrVMD 

gf X ) = e V^7 ^ 6 i_^P (2?) 

4i^ s - m 2 b + im^Y^il - x)s - m 2 + impY p ((l - x)s) 



with 371,2 = ^ 1 anc i x = 2 — x\ — X2- The quantities g^, g% are the coupling 

S 



constants determining the — > pit and p — > 7T7 verteces, respectively, 



'3T(0 -> e+e~) 



and e = v^a. A fit to the KLOE data — > 7r°7r°7 [2H] gives: 



^ = 0.811G6V- 1 , < 7 = 0.295GeV-\ = 42.5, nJ MD = 0.58195, (28) 

/3 P = 32.996°, (3^ = 163°. 



4 Including n^ A/£) one rescales the coupling constant. In our opinion it rescales the constant that cannot 
be directly determined from any experimental decay width. 



6 



2.4 Other contributions. 

We included in our program the channel 7* — » p ± n T — > 0, whose amplitude has been 

evaluated in RPT model. To write this part of FSR we used the results of Ref. |20j for the 
function f? : 



A/f 



8HI 



{k-Q + l 2 ) 
1 



1 



1 



C(l) C(-l) 



2k- 1 



1 



1 



C(l) C{-1) 



QAHy 

~9/r 



+ 



sh 2 



9/2 v<?(0 (7(-o 
1 1 



9/2 \C{1) C(-l) 



(29) 
(30) 
(31) 



where C{±1) 



m. 



(k + p ± ) 2 -tm p T p ((k+p ± ) 2 ) with (k + p ± ) 2 = (Q 2 + l 2 + 2k-Q±4k-l)/4. 



A value of the constant Hy is determined by the width of the p — > 717 decay 

r( p ± - vr± 7 ) 



Aam 3 p H v 



and can be related with the constant : if' 



27/2 



3f«g p 



1 



7T\3 



7T7 



. It gives Hy = 0.0144. In agreement 



with the calculation given in [20J, we found a negligible contribution of this channel and, for 
the sake of simplicity, we discard its effects on the numerical results presented in the following 
section. 



3 Numerical results 

In this chapter we present the results for the differential cross section and the forward-backward 
asymmetry [12], [29] for the reaction e + e~ — ► 7r + 7r _ 7 , where the FSR amplitude (Mfsr) receives 
contributions from both RPT (M R pt) and the <fr — ► 7r + 7r _ 7 decay (M^). The last one is a sum 
of the (f) direct decay (M^ cal ) and VMD (MY MD ) contributions. Thus the total contribution of 
the process ([2]) dar can be written as: 

da T = da! + da F + da IF ~ \M ISR + M F sr\ 2 , (32) 

dai ~ \M ISR \ 2 , 

da F ~ \M RPT \ 2 + \M (P \ 2 + 2Re{M R p T -M;}, 

da IF ~ 2Re{M /5iJ ■ {M RPT + M^)*} 



where daj corresponds to the ISR cross section, da F is for the FSR one. The interference term 
da if is equal to zero for symmetric cuts on the polar angle of the pions [7j. 

The different contributions to the FSR differential cross section daF, evaluated at s = 
m^, are shown in Fig.l, right, for the full angular range 0° < # 7 < 180°, 0° < 9 n < 180°. 
Good agreement between the results of the Monte Carlo simulation (points), with the analytic 

5 In the energy region s < this direct transition 7* — * p ± n Zf can be considered as the tail of the double 
resonance contribution of the p' meson decay: 7* — ► p' — > pn for s = m 2 . 
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Figure 1: {Left) The dependence of the branching ratio of the direct decay on the intermediate 
scalar states. {Right) Contribution to the FSR cross section dop/dQ 2 in the region 0° < 9 1 < 
180°, 0° < 9 n < 180° at s = m^. RPT is represented by circles, sQED by crosses, by triangles, 
while the difference between RPT and sQED is indicated by squares. 

prediction (solid line) is found. It can be noted that at low Q 2 the contribution from the 
direct decay (i.e. the term proportional to \M^ cal \ 2 in Eq. ( 1321) ) is quite large and, therefore, 
the additional contribution beyond sQED, can be revealed only in the case of destructive 
interference between the two amplitudes {Re{M RPT ■ M^) < 0). Published data from the 
KLOE experiment [2] are in favour of this assumption, which we will use in the following. 

To begin with, we consider the case s = m 2 ^. In Fig. 2 we show the values of the differential 
cross section dap/ da i and the forward-backward asymetry for the angular cuts of the KLOE 
large angle analysis [T6J, [30] : 

50° < 6 1 < 130°, 

50° < 9 W < 130° (33) 

for the Bremsstrahlung FS process in the framework of the sQED*VMD model and with the 
decay contributions (VMD and the direct decay), for a hard photon radiation with energies 
£ 7 > 20 MeV. 

Figure 3 shows the effects of RPT and direct decay terms to the total differential cross 
section and their contribution to FSR one for angular cuts (|33|) . Several distinctive features can 
be noted: (1) the peak at about 1 GeV 2 corresponds to the f intermediate state for the direct 
— > 7T7T7 amplitude; (2) the presence of RPT terms in the FSR is relevant at low Q 2 , where 
they give an additional contribution up to 40% to the ratio do 'rpt^i 'do s qed+^ (as shown in 
Fig. 3, left); (3) the destructive interference between the decay an the Bremsstrahlung FS 
process decreases the visible cross section in the whole Q 2 region and its dependence on FS 
Bremsstrahlung model at low Q 2 (see Fig. 2, left, down). Also we would like to draw attention 
to the VMD contribution. As we can see from Fig. 2, the VMD contribution almost does not 
change the value of the differential cross section (Fig. 2, left), but it changes essentially the 
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Figure 2: The ratio of the total cross section with respect to ISR one dar/dai (left) and the 
forward-backward asymmetry (right) as function of the invariant mass of the two pions, when 
the 4> contribution is taken into account and the Bremsstrahlung process is in the framework of 
the sQED*VMD model. The angular region is 50° < 6 1 < 130°, 50° < 6 n < 130° and s = ml. 



2.6 
2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 

















rSK ; 


= syi;» 








» FSR = 


= sQKD+9« !jj 








• FSR : 


= Ri'T+o,:*'^ 

>^ 




























L 


,„ L.n---"''*."' 

■hhheZm****/... 




■ "| ii'mim-nnTT'" 
: , , , , i , , , , I , , . , 


, , , . 


.... i ... . 


.... i .... i ... . 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



2.2 



0.8 0.9 1 

Q 2 (GeV 2 ) 



e 

f 1.8 

1.6 
1.4 
1.2 
1 













o FSR = 


sOFX 


1 














FSR = 


iQEI>+(|) 














• F$R = 


RPT-j-<|> 






































-. . , , 


, , , , 






, , , , 













0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 0.3 0.325 0.35 

Q 2 (GeV 2 ) 



2.2 
2 
1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 



*u 






• C 


°RPt + / dC 

.«! -/J**.. 


sQED 




■ +++ 


*4 




~*-t 






)ED 






f++ 




























#- 

- -♦- j 
* 
















• 
































". . . . i . . . . 


. , . , i , . , , 




. , , , 






, , , , 


, , , , 



0.1 0.1250.150.175 0.2 0.2250.250.275 0.3 0.3250.35 

Q 2 (GeV 2 ) 



Figure 3: The ratio dax/dai (left) as function of the invariant mass of the two pions for 
different models describing the Bremsstrahlung FS process (either RPT or sQED*VMD) and 
the ratio of FSR cross section in the framework of RPT, with respect to sQED, when the <ft 
direct decay contribution is (or not) taken into account (right), in the region 50° < # 7 < 130°, 
50° < 9 V < 130° at s = ml. 



9 



3 
2.5 

2 
1.5 

1 

0.5 

0. 











if... 


ESR: 
FSR: 


- sQKI) 




A 












= sQED+(|) 




I 








• 


FSR: 


= RPt+Oi 




Dtp 












"A"' 


























i , , . , 


, , , . 




, , , , 






, , , , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Q 2 (GeV 2 ) 



T3 



2.5 
2 

1.5 
1 

0.5 



0, 













o F5 


?R ..=.. 


-.Oil) 




: -♦- 


>- 








FSR = 


sQED+d) 






• • 








• F5 


>R = 


RPT- 


















i- ■ m 




























, , , , 








, , , , 






, , , , 





0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 0.3 0.325 0.35 

Q 2 (GeV 2 ) 



2 
1.9 
1.8 
1.7 
1.6 
1.5 
1.4 
1.3 h 
1.2 













• d 


:./J«... 


sQED 






T 










U RP1 




)ED 




























..... 









































































, , , , 












. , , , 






, , , , 



0.1 0.1250.150.175 0.2 0.2250.250.275 



0.3 0.3250.35 

Q 2 (GeV 2 ) 



Figure 4: The same ratio as in Fig. 3 at s = lGeV 2 . 



value of the forward-backward asymmetry (Fig. 2, right) and it follows the experimental data 
for it. Last but not least, all contributions beyond sQED are large enough near the threshold 
to make the analysis difficult. 

In order to reduce the background from <fi decay in the measurement of the pion form 
factor at threshold, KLOE has taken more than 200 pb _1 of data at 1 GeV [31]. In this case 
the meson intermediate contributions are suppressed (see Fig. 4, (right), the values of dar 
with and without the <j) decay almost coincide) and the main contribution to FSR comes from 
Bremsstrahlung process (see Section I2~TT1) . allowing to study a model for it. 

4 Model-independent test of FSR models 

Contributions to Bremsstrahlung FS process beyond sQED, as in the case of RPT, can lead to 
sizeable effects on the cross section and asymmetry at threshold, as shown in Figs. 2-4. Precise 
measurement of the pion form factor in this region needs to control them at the required level 
of accuracy. This looks like a rather difficult task, if one thinks that effects beyond sQED, as 
well as the contribution from — > 7r + 7r~7, are mo del- dependent. 

One can think to construct a general amplitude for the e + e~ — > 7r + 7r~7, according to some 
underlying theory, and try to determine the free parameters by a constrained fit on specific 
variables (like mass spectrum, charge and forward-backward asymmetry, angular distribution, 
etc.). Particularly for the charge asymmetry, it has been proved to be a powerful tool to 
discriminate between different models of (f) — > 7r + 7r~7 [29]. However, when the number of the 
parameters is large, correlations between the parameters of the model can arise and spoil the 
effective power of these fits. The situation becomes even worse if the pion form factor also has 
to be extracted from the same data. As an example, in the case of RPT model, if we consider 
only the p and u contribution to the pion form factor and the p and a\ contribution to FSR 
the number of free parameters is already six. The presence of the <fi direct and VMD decays 
results in additional free parameters. 
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The possibility to determine some of the parameters by external data can strongly help, as in 
the case of the — > 7r + 7r~7 amplitude, which can be determined by the 7r°7r°7 channel copiously 
produced at DA$NE. An additional source of information, which will be used to determine the 
contributions to FSR beyond sQED in a model-dependent way, is the dependence of the FSR 
amplitude on the e + e~ invariant mass squared s. 

Let us write the differential cross section for the emission of one photon in the process 
e + e~ — > 7r + 7r~7 as a function of the invariant mass of the two pions: 

$). = MW + ®, (34) 

where is the differential cross section for the emission of a photon in the final state, 

while the ISR function H S (Q 2 ) was defined in the Introduction. We indicate by a subscript s 
the dependence of each quantity on the e + e~ invariant mass (s). Since we will consider only 
symmetric angular cuts for pions, the interference term between initial and final state radiation 
has been neglected. 

At relatively high Q 2 the FSR differential cross section, , is dominated by the 

contribution coming from sQED (M s qed) and direct decay (M^): 

(^^) r \M sQED + M,\ 2 . (35) 

Contributions beyond sQED (AM) are expected to be important at low Q 2 . They introduce 
an additional term (AM) in the above expression: 



/ dc f \ 
\d~Q 2 ) ! 



\M sQED + AM + M^| 2 = (36) 
\M sQED + M/ + |AM| 2 + 2Re{AM • (M sQED + M^)*|. (37) 



We will now consider the following quantity: 

/ da T \ / df_sQE_D+±\ 



Y S (Q 2 ) = VQ2)s H \ Q ^ 2 )s = \FAQ 2 )\ 2 + AF S (Q' 2 ), (38) 

where AF S ~ (|AM| 2 + 2Re jAM • (M sQED + Mf)*}} /H 8 . 

If no contribution beyond sQED is present (AM = 0), Y S {Q 2 ) coincides with the square of 
the pion form factor, independently of the energy y/s at which it is evaluated, while any 
dependence on s is only due to an additional contribution to FSR beyond sQED. In particular, 
the difference of Y S {Q 2 ) computed at two beam energies {s\ and S2), can only come from FSR 
beyond sQED: 

AY(Q 2 ) = Y n (Q 2 ) - Y S2 (Q 2 ) = AF Sl (Q 2 ) - AF S2 (Q 2 ) (39) 

Therefore, before extracting the pion form factor at threshold, we suggest to look at the dif- 
ference AY(Q 2 ), which can be used to estimate the contribution beyond sQED to the FSR 
amplitude in a model-independent way. 

As a realistic application of this procedure, we consider the case of DA$NE, where KLOE 
has already collected more than 200 pb _1 at 1 GeV 2 and 2.5 fb" 1 at m^, that, in the range 
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Figure 5: Left: Y S (Q 2 ) at s = 1 GeV 2 (triangles), and at s = m| (circles), when FSR includes 
only the sQED and <fi contribution. The pion form factor \F n (Q 2 )\ 2 is shown by a solid line. 
Right: The difference AY(Q 2 ). 
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Figure 6: Left: Y S (Q 2 ) at s = 1 GeV 2 (triangles), and at s 



(circles), when FSR includes 



RPT and <fi contribution. The pion form factor \F n {Q 2 ) 
difference AY(Q 2 ). 
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is shown by a solid line. Right: The 
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Q 2 < 0.35 GeV 2 , correspond to O(10 3 ) and O(10 4 ) events, respectively, in the region 50° < 
# 7 < 130°, 50° < 6 V < 130°. We will consider RPT as a model for the effects beyond sQED. 

Fig. 5, left, shows the quantity Y S {Q 2 ) at S\ — 1 GeV 2 and at s 2 = m 2 when no additional 
RPT term is included in FSR. As expected, each of these quantities coincides with the square 
of the pion form factor \F n (Q 2 )\ 2 , shown by a solid line. The difference AY(Q 2 ) is shown in 
Fig. 5, right, which is consistent with zero as expected. In the region s < 0.35 GeV 2 one can 
expand the pion form factor as it was done in [I]: 

F n (q 2 ) ~1 +Pl . q 2 +P2 . q \ (40) 

Using the same experimental data for the pion form factor [1] as before we have: p\ = 1.15±0.06 
GeV~ 2 , p 2 = 9.06 ± 0.25 GeV" 4 , ^jv ~ 0.13. A combined fit of Y S (Q 2 ) to the pion form factor 
gives the following values: p 1 = 1.4 ± 0.186 GeV" 2 , p 2 = 8.8 ± 0.73 GeV" 4 , ^jv = 0.25, that 
is in a reasonable agreement with the results (140p . 

The situation is different as soon as the Bremsstrahlung FS process is modeled by RPT. In 
this CclSC, clS shown in Fig. 6, right, the difference AV(Q 2 ) ^ and the quantities Y S (Q 2 ) cannot 
be anymore identified with \Ft T (Q 2 )\ 2 , (see Fig. 6, left) U- 

Before concluding, we would like to point out the main points of our present method: 

• The quantity ( ~~ £ §j|ir ± ^' J is an input parameter of our procedure, and can be computed 
numerically by Monte Carlo; 

• The amplitude for <fi — > 7r + 7r _ 7 is taken from the 7r°7r°7 channel. 

• missing ISR multi-photon radiative correction can be added in H s and it will not spoil 
the effective power of the method; 

• A clear advantage of the procedure based on a Monte Carlo event generator is that it 
allows to keep control over efficiency and resolution of the detector and fine tuning of the 
parameters. 

Even if the main limitation of the method could come from the uncertainty on the parameters 
of — > 7r + 7T - 7 amplitude, especially at low Q 2 , we believe that the KLOE data on <fi — > 7r°7r°7 
will allow a precise description of this amplitude. In any case, in agreement with |29j we strongly 
recommend to check the amplitude by using charge asymmetry and to compare with spectrum 
of the 7r + 7r~7, at least at high Q 2 , where the pointlike approximation is safe (as done in |30j). 

5 Conclusion 

A test of FSR at threshold in the process e + e~ — > 7r + 7r~7 is a rather important issue, not 
only for the role of FSR as background to the measurement of the pion form factor, but also 
to get information about pion-photon interactions when the intermediate hadrons are far off 
shell. At s — m 2 an additional complication arises: the presence of the decay <p — ► 7r + 7r _ 7, 
that goes either through the intermediate scalar (the direct decay) or the vector state (VMD 
contribution) whose amplitude and relative phase can be described according to some model. 

destructive interference between RPT and <f> — > ir + ir~j amplitudes tends to cancel out the effects beyond 
sQED at s = rr& (see Fig. 3, left). Therefore the quantity Y S (Q 2 ) almost coincides with the pion form factor. 



13 



By means of a Monte Carlo event generator FEVA, which also includes the contribution of the 
rare <p ~> 7r + 7r~7 decay, we estimate the effects beyond sQED in the framework of Resonance 
Perturbation Theory (RPT) for angular cuts used in the KLOE analysis of the pion form factor 
at threshold. We show that the low Q 2 region is sensitive both to the inclusion of additional 
terms in the FSR amplitude given by the RPT model and to the <fi decay contribution (especially 
its VMD part). 

We also propose a method which allows to estimate the effects beyond sQED in a model- 
independent way. We found that the deviation from sQED predicted by RPT can be observed 
within the current KLOE statistics. 

We would like to emphasize once again that this work was motivated by ongoing experiment 
on precise measurements of the muon anomalous magnetic moment [32] that allows to perform 
tests of the Standard Model with a fabolous precision. 
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